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Abstract 

We calculate the thermodynamic potential of bubbles of hadrons embedded 
in quark-gluon plasma, and of droplets of quark-gluon plasma embedded in 
hadron phase. This is a generalization of our previous results to the case of 
non-zero chemical potentials. As in the zero chemical potential case, we find 
that a quark-gluon plasma in thermodynamic equilibrium may contain stable 
bubbles of hadrons of radius R ~ 1 fm. The calculations are performed within 
the MIT Bag model, using an improved multiple reflection expansion. The re- 
sults are of relevance for neutron star phenomenology and for ultrarelativistic 
heavy ion collisions. 
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I. INTRODUCTION 



In a previous paper [fl] we calculated the free energy of a bubble of hadrons embedded in 
an extended quark-gluon plasma (QGP), and of a droplet of QGP embedded in an extended 
hadron phase, for parameters in the vicinity of the cosmological quark-hadron transition, i.e. 
the baryon chemical potential were set to zero. In Ref. |J we found, as Mardor and Svetitsky 
in HH, that the free energy of a hadron bubble of radius R embedded in QGP possessed 
a minimum at radii of a few fm, even above the phase transition temperature. Thus, 
hadronization from QGP is strongly enhanced compared to the usual nucleation scenario, 
where an energy barrier has to be overcome before bubbles of hadrons can form and grow. 

An important ingredient in these calculations is the density of states of the relevant par- 
ticles. In Ref. H we advised a modification of the usual expressions arising from the multiple 
reflection expansion (MRE), and we saw that this modification (the MMRE) yielded more 
accurately the free energy than the MRE did when compared to a direct (but numerically 
heavier) sum-over-discrete-states calculation. Our modification of the multiple reflection 
expansion expressions is physically well motivated, in that it consists of a truncation of the 
density of states in such a way that we avoid the effects of a negative density of states, 
present in the usual MRE expressions. 

In this paper, we generalize the calculations to finite chemical potential, and we shall 
see that also in the case of non-zero chemical potential, the MMRE produces accurately 
the thermodynamic potential. We show that, as in the case of zero chemical potential, the 
thermodynamic potential of a hadron bubble embedded in QGP has a minimum at a radius 
of R ~ 1 fm, meaning that it is thermodynamically favorable for a QGP in equilibrium 
to spontaneously create bubbles of hadrons of about this size. However, the minimum in 
the thermodynamic potential found for non-zero chemical potential is less pronounced than 
in the zero chemical potential case considered in Ref. 0. The more general calculations 
presented in this paper should be relevant also to neutron star phenomenology as well as 
ultrarelativistic heavy ion collisions. 

In the following section we give an overview of the basic theory leading to the results of 
section fTJ. Finally in Sec. (TV| , we summarize the conclusions. 

II. THEORETICAL FRAMEWORK 
A. The particle model 

1. The quark-gluon plasma 

The model of the plasma phase is a close derivative of the MIT Bag model 0||. It 
consists essentially of non-interacting quarks and gluons (also self-interaction of the gluons 
is neglected), bounded by MIT Bag-like boundary conditions ||. There are three quark 
flavors in our model: u and d, which we consider massless, and s, to which we assign a mass 
of m s = 150 MeV. The other quarks are too heavy to be relevant in this analysis. We choose 
a critical temperature (the bulk phase equilibrium temperature at zero chemical potential) 
of T = 150 MeV, which gives a Bag constant of B = (221 MeV) 4 . The plasma may be inside 
a sphere, as in the usual MIT Bag (we shall refer to this configuration as a plasma droplet), 
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or outside, i.e. in a "vacuum bubble" configuration. For a detailed review of this model, the 
reader may consult We are mostly interested in the vacuum bubble situation, where the 
vacuum bubble is filled with a thermal mixture of hadrons as described later. 

For the density of states of the quarks and gluons we use a modification of the multiple 
reflection expansion (MRE) |6|-|T0||, advised in ]![. This modification affects only the gluon 
density of states, which is then 

/ 0, 0<kR< 0.832 

P*( k > R ) = [Yag*_u kR > . 8 3 2 C 1 ) 

for gluons inside a sphere (i.e. the MMRE(i2)), and 

{ 0, 0<kR< 0.458 

+ M ^> . 458 (2) 

for the difference between the gluon density of states in a volume and in volume — Vh, 
(i.e. the MMRE(-i?)). 

In |l[ we saw that this modified multiple reflection expansion (MMRE) description of 
the density of states made the free energy agree nicely with more direct calculations, i.e. 
summing over discrete energy levels. Before using the MMRE in the more general case 
under consideration in this paper, we performed similar checks of the MMRE against direct 
sum-over-discrete-states calculations. Also in this case, the MMRE proved to describe the 
density of states adequately and significantly better than the MRE, cf. Fig. [I]. For a complete 
description of the MMRE density of states, we refer the reader to 

2. The hadron phase 

We include all hadrons with masses below 1.2 GeV, taking only the volume part of the 
density of states into account. This is justified as the hadrons represent far fewer degrees of 
freedom than the QGP. Specifically, the density of states of a hadron species occupying a 
volume Vh = ^-R 3 is 

PH (k,R)=y^. (3) 



B. Thermodynamics 

The thermodynamic potential is 

tt(T, V, {/i,}) = ~T HZ(T, V, {/!,})), (4) 
where Z(T, V, {/ij}) is the partition function 

Z(T, V, {in}) = Tr { e -(w-£iMi)/Tj . (5) 
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Here, TC is the Hamiltonian, the Aj's are the symmetries of the Hamiltonian (the conserved 
quantities), and {/ij} denotes the corresponding collection of chemical potentials. Tr{..} 
means the sum over all diagonal elements of energy- and number-eigenstates. These are 
simultaneous eigenstates since [TC, Aj] = 0. The volume dependence enters via the energy 
levels, i.e. in the trace operation. 

In thermodynamic equilibrium, the configuration realized in Nature is the one which min- 
imizes the thermodynamic potential, leading to the Gibbs conditions for phase equilibrium: 
T\ = T 2 (thermal equilibrium), — |^ = —^r (mechanical equilibrium), and /4 = /if 1 ' 
(chemical equilibrium), where index 1 and 2 refer to the two phases. 

If there are no interactions between the particles, i.e. the energy of a particular state is 
the sum of the energies of the individual particles, then we can write the thermodynamic 
potential for a particular particle species enclosed in a spherical volume of radius R, as 

n(R,T,fj,) = T9T [°° dkp(k,R)\n(l±e~ iVWT ^-^ /T ) 7 (6) 
Jo 

(upper sign for fermions, lower sign for bosons) where we further assume that the energy 
levels are sufficiently closely spaced that we may use a smoothed density of states, p(k,R), 
(e.g. the MMRE) normalized such that J °° dk p(k, R) is the total number of states in the 
volume V = ^-R?. The chemical potential entering in (||), is the combined chemical potential 
of that particle species 

li = Kft, (7) 

i 

Aj being the quantum expectation value of the symmetry operator Aj of the particle species 
in question (e.g. Abaryon = 1/3 for a quark). Finally, g is an appropriate degeneracy factor. 

At the energies relevant in this analysis, the quarks and gluons can be considered point- 
like. This is not the case for the hadrons. We must take explicit account of the fact that the 
hadrons occupy a certain volume. A number of such excluded volume corrections have been 
discussed in the literature, but we shall not go into an analysis of the different suggestions. 



We choose the one proposed in []TT[ , since it is very easy to implement and grasps much of 
the essential physics involved. According to [fLlp , the "true" pressure of the hadron gas (pn) 
is obtained from the ideal gas pressure (pH,id) according to 

PH,id , Q s 

where 

i 

is the total energy density of the ideal hadron gas, 

PH,id = ^2PiM (!0) 

i 

is the total hadronic ideal gas pressure, (i G hadron species), and B is the Bag constant. We 
note that some sort of excluded volume correction is essential, since in a point-like hadron 
gas, no transition to QGP will occur at low temperatures, even at arbitrarily high density. 
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There is one more question about the hadrons we need to address before we proceed. 
This is the phenomenon of Bose-Einstein condensation, regarding the bosonic part of the 
hadron gas. A bosonic gas always has /i, < and when /Zj = m ; particles will start to 
"condense" into the lowest, zero-momentum state, meaning that there will be a macroscopic 
number of particles in this state. When the parameters are such that we would have fii > rrii, 
the hadron species i is removed from the above summations (|j) and (|T0|) , since particles in 
a zero-momentum state do not contribute any pressure. 

III. RESULTS 

We shall consider both the case of equal chemical potentials, fi u = fid = fi s = fi g , and 
the case where the u and d quarks have a different chemical potential from that of the s 
quark, /i u = /i^ 7^ /i s . The first case is relevant when the time scales are such that the weak 
interactions maintain equilibrium between all three flavors of quarks (e.g. quark matter in 
neutron stars, neglecting the electron chemical potential as a first approximation), whereas 
the case of separate s quark chemical potential is relevant to ultrarelativistic heavy ion 
collisions, the time scales here being such that the net number of s quarks is conserved 
separately from that of the net number of light quarks. In our model, the two light quarks 
are both massless, and since we neglect electromagnetic interactions, there is no difference 
between u and d quarks. 

A. Quark-hadron phase equilibrium 

To study the phase equilibrium between the hadrons and the QGP, we implement the 
Gibbs conditions on the volume (or bulk) part of the thermodynamic potential, i.e. also in 
the QGP density of states we only use the terms proportional to the volume. (We shall see 
later that the surface terms give rise to interesting features near bulk phase equilibrium.) 
In this way, we obtain the phase diagram in Fig. |^. The QGP phase (above the phase 
equilibrium lines) occupies a larger and larger region of the phase diagram as the s-quark 
chemical potential is raised. As one might have expected, the phase equilibrium line in the 
case of equal chemical potentials lies somewhere between the lines of varying fi s . 

B. The thermodynamic potential 

In this section we look at the thermodynamic potential in the vicinity of equilibrium, 
now using the full expression of the density of states of the QGP (i.e. including the surface 
terms). In Figs. [3]-|6|, there is a single chemical potential, fi q , common to all three quark 
flavors. 

First we show Fig. |3], the thermodynamic potential of a QGP droplet in equilibrium 
with an extended hadron phase near the bulk equilibrium point (fi q ,T) = (320, 106.1)MeV. 
The picture is as expected for a first order phase transition: Below the bulk equilibrium 
temperature no stable QGP droplet can form, and even somewhat above this temperature 
there is an energy barrier for the system to pass before the true minimum at R = 00 can be 
reached. The surface terms are responsible for this energy barrier. Precisely at the phase 
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transition temperature the volume terms cancel (by definition), and the thermodynamic 
potential goes to infinity as R? (the leading surface term). 

The more interesting conclusions are reached when one considers the reverse situation, 
namely a hadron bubble in equilibrium with an extended phase of QGP. Figs. |J and |6| 
show the thermodynamic potential of this configuration near different points on the phase 
equilibrium line. In all three cases the potential exhibits a minimum at a radius of R ~ 1 fm 
of approximately the same depth, even at temperatures above phase equilibrium tempera- 
ture. Such a minimum would not be present in the standard textbook treatment of phase 
transitions, in which a first order phase transition usually is described in terms of a phe- 
nomenological free energy containing only volume- and (positive) surface terms. 

Thus, a QGP in thermodynamic equilibrium apparently contains bubbles of hadrons, and 
the transition from QGP to hadrons will proceed (in addition to ordinary bubble nucleation 
and subsequent expansion of these bubbles) by smooth expansion of the pre-formed bubbles. 
In addition to this smooth growth as the minimum tends to larger radii when the temperature 
drops, pre-formed bubbles can also grow discontinuously by passing the barrier from R > 
to some R outside the barrier. This possibility resembles very much the ordinary nucleation 
scenario with a modified barrier height. At the parameters of Figs. |], |5] and || the mean 
number of quarks in a hadron bubble of radius R = 1 fm is typically 3-4, suggesting that 
formation of a real hadron is not a very rare event. However, the minimum being of modest 
depth, this interesting effect could be due to the model's inaccurate representation of QCD. 

The fact that there is no energy barrier for a hadron bubble to form in QGP does 
not imply that the phase transition is second order, or a smooth cross-over. In the model 
adopted here, the transition is inherently first order, since there is an a priori difference in 
entropy between the two bulk phases, and thus a latent heat. As seen from e.g. Fig. [| it 
is not energetically favorable for bulk hadronic matter to form at temperatures above the 
phase transition temperature To. But the formation of isolated bubbles of hadronic matter 
is favorable above To, so the transition will appear smoother than normally expected for a 
first order transition, with a more gradual release of latent heat. 

We now investigate the effect of letting the s quark have its own chemical potential. To 
this end we show Figs. [7] and |8|. In Fig. [7| the s quark chemical potential is 4 times greater 
than the light quark chemical potential, whereas in Fig. || it is converse. The conclusion to 
be drawn from these figures is that a large s quark chemical potential tends to wash out the 
minimum in the thermodynamic potential, present when the light quark chemical potential 
is less than or equal to the s quark chemical potential. This behavior can be traced back to 
the positive surface tension contribution from the massive s-quarks to the thermodynamic 
potential. 

Finally, we underline that the interesting effect of bubble formation is due to the fact 
that the presence of a surface alters the density of states of particles. Although these finite- 
size effects are negligible at large radii of the bubble or droplet, they may have important 
implications for the phase transition as a whole. 

IV. CONCLUSION 

We have seen that the thermodynamic potential of a hadron bubble embedded in quark- 
gluon plasma, exhibits a minimum at a radius R ~ 1 fm, even at temperatures somewhat 
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above the bulk transition temperature. Thus, within the model described here, a homoge- 
neous plasma of size larger than a few fermi is an impossibility. 

Regarding the relevance of these results in connection with current and forthcoming 
ultrarelativistic heavy ion collisions, we conclude that if indeed a quark-gluon plasma is 
formed in the course of such a collision, then (according to the model considered here) this 
plasma phase will contain bubbles of hadrons of radius R ~ 1 fm. If such hadronic bubbles 
form inside the plasma, the observable effects are likely to include some blurring of the 
plasma signatures, since this reduces the effective plasma volume. 

We emphasize that these conclusions may well be model dependent, inasmuch the min- 
imum of the thermodynamic potential is of modest depth. However, the phenomenon of 
spontaneous creation of stable hadronic bubbles in a quark-gluon plasma, does seem to be 
well established within the the model discussed here fl|,3, |1^1 ■ The main question is now, 
whether these hadronic bubbles are of physical nature, or merely an artifact of the model. 
Certainly, it would improve confidence in these results if other models of QCD, and eventu- 
ally lattice calculations, were to yield similar results. 
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FIG. 1. Thermodynamic potential of a bubble of hadron phase of radius R embedded in an 
extended QGP phase. In the two upper curves (fi q ,T) = (450,46) MeV, i.e. a higher temperature 
than the phase transition temperature at this chemical potential, To = 41.0 MeV. In the two 
lower curves (fJ. q ,T) = (200,139) MeV, again slightly above the phase transition temperature 
To = 133.7 MeV. To draw the solid lines, the MMRE(— R) has been used to describe the density 
of states in the expression for the thermodynamic potential, whereas the dashed lines are for the 
MRE(— R). Based on comparisons with more direct sum-over-states like calculations (the points 
practically coinciding with the MMRE(— R) lines), we conclude that the MMRE(— R) is the more 
correct model for the density of states. 
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FIG. 2. Phase equilibrium lines in the (fj, u , T)-plane. At these lines, QGP and hadrons are in 
thermodynamic equilibrium. The solid line represents the case of equal quark chemical potentials, 
fJ'u = t^d = Ms> whereas the other lines represent different values of the s-quark chemical potential: 
Top to bottom, fi s = 0, 100, 200, 300, 400 MeV. 
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FIG. 3. Thermodynamic potential of a droplet of QGP of radius R embedded in an extended 
phase of hadrons, normalized such that the thermodynamic potential of a pure hadron phase is 
zero. The plasma density of states is described by the MMRE(i?). The bulk phase transition point 
is (fj, q ,T) = (320,106.1) MeV (cf. Fig. |2|). Curves are shown for different temperatures, top to 
bottom: T = 104, 106.1, 108, 112 MeV. The chemical potential is fixed at fi q = 320 MeV. 
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FIG. 4. Thermodynamic potential of a bubble of hadron phase of radius R embedded in an 
extended QGP phase, normalized such that the thermodynamic potential of a pure QGP phase is 
zero. The plasma density of states is described by the MMRE(— R). The bulk phase transition 
point is (jiq,T) = (320, 106.1) MeV (cf. Fig. |). The chemical potential is fixed at fi q = 320 MeV, 
and the temperature varies, top to bottom: T = 112, 108, 106.1, 104 MeV. 
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FIG. 5. As Fig. ^, but this time investigating the bulk phase transition point 
(fj, q ,T) = (450,41.0) MeV. The chemical potential is fixed at this value, and the temperatures 
are, top to bottom: T = 52,44,41.0,38 MeV. 
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FIG. 6. As Fig. ^, but this time investigating the bulk phase transition point 
(Hq,T) = (470.0,3.28) MeV. The chemical potential is fixed at this value, and the temperatures 
are, top to bottom: T = 8, 3.28, 1 MeV. Notice that the depth of the minimum at R ~ 1 fm is 
practically independent of which equilibrium point we consider, cf. Figs. || and|B[ The curves in this 
figure differ somewhat in their qualitative behaviour from those in Figs. || and || This is because 
at these low temperatures, the surface terms dominate over the volume terms at the radii shown 
here; at larger radii the T = 1 MeV curve bends over and eventually tends to — oo as it should. 
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FIG. 7. Thermodynamic potential of a hadron bubble embedded in an extended QGP phase 
in the vicinity of the phase equilibrium point \i u = 100 MeV, fi s = 400 MeV, To = 119.5 MeV 
(cf. Fig. [2]). The chemical potentials are fixed at these values, and the temperature varies, top 
to bottom: T = 125, 119.5, 115 MeV. Notice that, at these values of the chemical potentials, the 
minimum in the thermodynamic potential practically disappears. 
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FIG. 8. As Fig. 0, but the chemical potentials are now fj, u = 400 MeV fi s = 100 MeV, and 
the phase equilibrium temperature corresponding to these values is To = 100.0 MeV. Again the 
chemical potentials are fixed, and the temperatures are, top to bottom: T = 105, 100, 90 MeV. The 
minimum in the thermodynamic potential is quite deep at these values of the parameters. 



12 



